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Abstract 


\c 


A one-dimensional  model  consisting  of  conserva- 
tion equations  for  momentum,  moisture,  and  thermody- 
namic energy  in  the  cloud-free  marine  atmosphere  is 
constructed  and  used  to  investigate  the  effect  of 
sea  surface  temperature  perturbations  on  the  local 
thermodynamic  characteristics  of  the  atmosphere. 

It  is  found  that  for  perturbations  in  the  range  that 
would  be  produced  by  Ocean' Thermal  Power  Plants,  the 
average  temperature  and  specific  humidity  in  the 
marine  boundary  layer  vary  linearly  with  the  sea 
surface  temperature  with  slopes  of  0.8  'C/?C  and 
0.6  £/icg  / C respectively.  The  total  heat  flux  at 
the  sea  surface  is  also  found  to  vary  almost  linear- 
ly with  the  sea  surface  temperatur^in  the  range 
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with  a slope  of  31  cal/c 
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Section  1 
INTRODUCTION 

Adequate  planning  for  the  deployment  of  large  numbers  of 
ocean  thermal  power  plants  requires  knowledge  of  their  effect 
on  the  thermodynamic  characteristics  of  the  marine  atmosphere. 
OTPP's  will  alter  these  characteristics  because  their  opera- 
tion will  necessarily  produce  a lowering  of  the  sea  surface 
temperature. 

It  is  essential  to  estimate  the  magnitude  of  this  environ- 
mental impact  for  several  reasons.  For  example,  the  rate 
at  which  the  thermal  resource  is  regenerated  downstream  from 
an  OTPP  depends,  in  part,  on  the  total  heat  balance  at  the 
sea  surface  which  is  a function  of  both  the  sea  surface  tem- 
perature and  the  thermodynamic  structure  of  the  atmosphere 
above.  An  estimate  of  this  regeneration  rate  is  necessary 
to  insure  that  OTPP's  are  not  positioned  where  they  will 
adversely  affect  other  plants  downstream.  Furthermore,  many 
OTPP's  will  probably  be  located  near  shore.  Hence,  changes 
they  ultimately  induce  in  the  thermodynamic  structure  of 
the  atmosphere  may  affect  the  climate  of  the  coastal  zone  by 
altering  the  dynamic  and  thermodynamic  characteristics  of  the 
sea  breeze  circulation.  Finally,  the  results  of  Shukla  (1975) 
suggest  that  widespread  sea  surface  temperature  anomalies 
may  significantly  affect  the  climate  of  regions  far  removed 
from  the  anomalies  by  reducing  evaporation  and  altering  the 
large-scale  atmospheric  flow  pattern.  Shukla  performed 
experiments  with  a general  circulation  model  and  found  that 
an  average  sea  surface  temperature  anomaly  of  -1.5°C  in  the 
Arabian  Sea  over  a region  roughly  the  size  of  the  Gulf  of 
Mexico  reduced  the  mean  rainfall  rate  in  India  by  about  50%. 
Such  results  are  disturbing  since  the  Gulf  of  Mexico  (a 
proposed  site  for  deployment  of  large  numbers  of  OTPP's)  is 
also  the  source  region  for  most  of  the  rainfall  in  the  mid- 
western  United  States. 
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Investigation  of  changes  induced  in  the  three-dimensional, 

r 

time- depen dent  state  of  the  atmosphere  requires  a three- 
dimensional,  time-dependent  model.  However,  studies  to  assess 
the  effect  of  OTPP  operation  on  the  local  time-averaged  state 
of  the  atmosphere  can  be  accomplished  with  a one-dimensional, 

**  steady  model.  Such  a model  has  been  constructed  for  the 

special  case  of  cloud-free  conditions  and  these  studies  are 
in  progress.  The  remainder  of  this  report  will  describe  the 
model,  some  preliminary  results,  and  plans  for  future  work. 


DESCRIPTION  OF  THE  MODEL 


BASIC  MODEL  EQUATIONS 


Consider  a region  above  the  tropical  ocean  with  an  area 
large  enough  to  encompass  the  sea  surface  temperature  anomaly 
produced  by  several  OTPP's  yet  small  enough  to  be  regarded  as 
a small  fraction  of  the  atmosphere's  general  circulation 
system.  Then,  in  the  absence  of  clouds,  the  equations  for 
conservation  of  momentum,  moisture,  and  thermodynamic  energy 
averaged  horizontally  in  the  region  can  be  written 


where 


has  been  used  as  the  vertical  coordinate.  All  the  symbols  .are 
defined  in  the  Appendix  and  most  of  the  notation  is  standard. 
The  first  terms  on  the  right-hand  sides  of  (1)  - (4)  represent 


the  effect  of  advection  by  the  average  vertical  motion  in 
the  region  while  the  second  terms  on  the  right  of  (Z  and 
(4)  represent  the  effect  of  horizontal  advection  of  moisture 
and  temperature  respectively.  The  treatment  of  the  hori- 
zontal moisture  and  temperature  advection  will  be  discussed 
.in  Section  2.4.  Note  that  the  Rossby  number  of  the  flow  is 
assumed  «1,  which  allows  neglect  of  the  horizontal  advection 
terms  in  equations  (1)  and  (2).  The  second  terms  on  the 
right-hand  sides  of  equations  (1)  and  (2)  represent  the 
Coriolis  force  and  the  horizontal  pressure  gradient.  The 
third  term  on  the  right  of  equation  (4)  represents  compres- 
sional  heating  while  the  fourth  term  is  the  net  radiational 
heating.  The  last  terms  on  the  right-hand  sides  of  equations 
(1)  to  (4)  represent  the  divergence  of  the  vertical  turbulent 
fluxes  of  E-W  momentum,  N-S  momentum,  moisture,  and  heat 
respectively. 


2.2  TURBULENCE  PARAMETERIZATION 

Following  many  others,  we  divide  the  planetary 
boundary  layer  into  two  parts: 

• a surface  or  "constant  flux"  layer,  with 
a thickness  of  order  20m,  in  which  the 
turbulent  fluxes  of  momentum,  moisture, 
and  heat  depart  little  from  their  surface 
values,  and 

• an  overlying  Ekman  layer  in  which  the 
fluxes  generally  decrease  with  height. 
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2.2.1  The  Surface  Layer 

Under  steady,  horizontally  homogeneous  conditions, 

Mon in-Obukhov  similarity  theory  predicts  that  properly 
scaled  vertical  derivatives  of  wind  speed,  moisture,  and 
temperature  will  be  universal  functions  of  z^/L  in  the  sur- 
face layer,  where  z ^ is  the  height  above  the  surface  and  L 
is  the  Monin-Obukhov  length.  We  impose  the  value  z1  « 20  m 
(the  height  of  the  first  mesh-point  in  our  representation.) 
The  length  L is  negative  in  unstable  cases,  as  described 
below.  The  surface  eddy  fluxes  are  then  related  to  the  wind 
speed,  temperature  and  humidity  at  z^ , according  to  the 
"profile  method"  (Busch,  1977),  by  the  following  equations: 


(v».>*  - «.  - *<u  - V 

!n(Zi/zo)  - \Pm(z1/L) 


H 


cp<u*(T  - Tq) 
lnUl/z0)  - iMz.,/L) 


6'~1' 


- qQ) 

ln(z1/zQ)  - H/q(z1/L) 


(5) 

(6) 

(7) 


Here 


Zj/L 

♦a  ■/  (i  - V‘»T 


for  a * m,0,q,  (8) 


where  $m,  4>0,  and  41^  are  non-dimensional  gradients  of  wind, 
speed,  temperature,  and  humidity  respectively.  Equations 
(6)  - (8)  have  been  used  recently  by  Hsu  (1974),  Pond,  et.  al. 
(1974),  Dunckel,  et.  al.  (1974),  and  several  others.  The 
forms  of  the  $a  functions  have  recently  been  reviewed  by  ' 
Hogstrom  (1974)  and  it  is  generally  accepted  that  most 
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atmospheric  data  are 

well  represented  by 

V')  - jo  - 
0 

150  * 

for  OO 
for  C<0 

(9) 

♦9  • ♦,  - Li ! 

90  9 

for  UO 
for  C<0 

(10) 

which  are  the  forms  used  in  the  model.  The  Monin-Obukhov 
length  is  calculated  from 


L 


pau* 


Kg(HQ/cpT  + 0.6Qq) 


(11) 


and  is  negative  when  the  surface  layer  is  hydrostatically 
unstable  (see  Kraus,  1972).  Following  Clarke  (1970)  and 
others,  we  assume 


zQ  = max  (0.032  u*2/g  , zmin)  , (12) 

where  z . is  0.0015  cm. 
mm 

The  sea  surface  temperature  provides  the  lower 
bou'dary  condition  for  T and  the  specific  humidity  at  the 
lower  boundary  is  set  equal  to  0.98  times  the  saturation  value 
at  the  sea  surface  temperature  (to  account  for  the  effect  of 
salt  on  the  vapor  pressure  of  water,  see  Roll,  1971).  The 
wind  velocity  at  the  sea  surface  is  set  equal  to  the  surface 
drift  velocity,  which  is  assumed  to  be  in  the  direction  of 
the  geostrophic  wind  at  the  top  of  the  Ekman  layer  and  have 
a magnitude  of 

°o  - ug  (P./0.X  w» 

-V  Ug/200 
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(see  Kraus,  1972).  Equations  (5)  - (12)  are  solved  iteratively 
at  each  time  step  in  the  numerical  representation  which  we 
have  adopted. 


2.2.2  The  Ekman  Layer 

The  vertical  eddy  fluxes  of  momentum,  heat,  and 
moisture,  in  the  Ekman  layer  are  calculated  from 
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and  the  virtual  temperature  Tv  is 


Ty  - T (1  + 0 . 61q ) 


7 


Following  Clarke  (1970),  the  mixing  length,  Si,  is  specified 
as  a function  of  z in  a manner  that  depends  on  the  sign  of  L. 
For  negative  L we  define  a "convection  depth"  height  z c£j 
above  which  the  Richardson  number 


Ri  = N2/)3u/3z)2  + ( 3v/3z)2 


is  greater  than  unity, 
frequency . 


Here  N is  the  usual  Brunt-Vaisala 


Then 


l = KZ(1  - 15  z/L)0,275  , for  z<0.3z 


l - l = 30  m 

o 


, for  z>z 


and  £ is  obtained  by  linear  interpolation  in  the  intervening 
range.  For  positive  L,  we  define  the  mixing  length 


£ = min  |<z(l  + 5z/L)_1  , £Qj  , for  z<L  , 
£ = min  jtcz/6  , £Q  j > f°r  > 


where  £q  is  again  30  m. 

This  treatment  of  the  vertical  diffusion  coefficient 
has  been  used  recently  by  Physick  (1976)  in  a two-dimensional 
model  of  the  sea  breeze. 
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2.3  RADIATIVE  TRANSFER  EQUATIONS 

The  long-wave  radiation  equations  used  in  our  model 
are  basically  due  to  Danard  (1969).  For  a cloud-free  atmos- 
phere, the  zero  approximation  transparency  yields  an  upward 

4 

energy  flux  of  eQB0T  . The  infrared  emissivity  e0  of  the 
sea  surface  is  taken  as  0.97.  The  first  approximation  to 
the  net  downward  long-wave  radiation  flux  at  a particular 
level  is  then 

* 4 C 4 

Ft  = - e BT  + / BT  de 
L 00  J 

above 

+ J U0BTo4  - BT4  )de . 

below 

Here  de  is  the  emissivity  and  absorptivity  of  the  layer 
between  the  level  where  the  temperature  is  T and  the  level 
in  question.  The  layers  above  add  downward  radiation 
according  to  their  temperature,  while  the  layers  below 
absorb  some  of  the  upward  radiation,  and  also  emit  new 
upward  radiation  according  to  their  temperature.  By 
definition,  de  is  given  by 

de  = e'(w)  dw 


where 


dw  = q (p/pr)°'85(Tr/T)0,5(p*/g)  do 

is  the  "precipitable  water  vapor"  in  the  layer  da  and 
c'  is  the  proportionality  derivative.  We  approximate 
c'  by  adapting  Elsasser’s  (1960)  experimental  curve  for 
the  flux  emissivity  of  pure  water  vapor. 
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The  incoming  short-wave  radiation  flux  is 
S^cos  a,  where  is  the  solar  constant  and  a is  the 
solar  zenith  angle.  This  flux  is  divided  into  two  parts, 
following  Arakawa,  Katayama  and  Mintz  (1968).  The  first 
part  (blue,  with  wavelength  smaller  than  about  0.9  y) 
experiences  significant  Rayleigh  scattering,  but  no 
absorption,  while  the  second  part  (red,  with  larger  wave- 
lengths) experiences  significant  absorption  but  no 
scattering.  The  net  downward  short-wave  radiation  flux 
distribution  at  a particular  level  is  taken  in  this  model 


651(1  - ao)/(l  - aoag) 


F’  - S0  cos  «[o. 651(1  - »0)/(l  - 

* 0 . 349 ] 1 - 0 . 03(w*  sec  a)0'303)] 


where  w is  the  integral  of  the  precipitable  water  vapor 
from  the  level  in  question  to  the  top  of  the  atmosphere, 
a is  the  albedo  of  the  sea  surface,  taken  as  0.065,  and 

S 


aQ  * 0.085  - 0.245  log 


I — cos  a I 
\PT  / 


where  pT  is  the  pressure  at  the  top  of  the  model. 

The  solar  zenith  angle  a is  a function  of  the  time 
of  day,  for  each  latitude  4>  and  solar  declination  angle  6 
from  the  equatorial  plane.  The  above  expression  for  the 
downward  flux  is  therefore  Integrated  over  the  daylight 
hours  to  obtain  the  diurnal  average  flux  distribution,  as 
a function  of  $ and  5. 
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2.4 


TREATMENT  OF  HORIZONTAL  MOISTURE  AND  TEMPERATURE 
ADVECTION 


The  horizontal  advection  of  moisture  averaged 
across  the  model  domain  is  parameterized  by 


JQ  _ 1)  (q  - qus> 


where  qug  is 


a specified  function  of  altitude  and  is  con- 
sidered to  be  the  specific  humidity  at  some  distance  of 
order  d upstream  from  the  center  of  the  region  being 
modeled  and  U is  the  wind  speed  at  the  level  in  question. 
Hence,  with  respect  to  moisture  advection,  the  model  can 
be  thought  of  as  representing  one  horizontal  grid  point  of 
a general  circulation  model  with  the  vertical  moisture 
profile  held  constant  at  the  grid  points  upstream  and  the 
lateral  advection  of  moisture  computed  with  an  upstream 
difference.  In  the  context  of  the  present  study,  qus  can 
be  identified  with  the  specific  humidity  distribution  just 
upstream  of  the  OTPP-generated  sea  surface  temperature 
anomaly. 

The  horizontal  advection  of  temperature  averaged 
across  the  model  domain  follows  from  the  thermal  wind 
equation  in  o coordinates  which  relates  the  horizontal 
gradient  of  temperature  to  the  geostrophic  wind  (which  is 
an  imposed  quantity  in  our  model).  Thus,  for  the  case  of 
flat  surface  topography,  we  have 


u ax  v ay  r 1 u i 

» - ff  n]| 


3T 

To 


with  all  the  symbols  defined  in  the  Appendix. 
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2.5  GRID  AND  FINITE  DIFFERENCE  SCHEME 

In  order  to  resolve  the  surface  layer  without 
demanding  unnecessary  core  storage,  a vertically  stretched 
grid  is  employed  with  the  stretching  function  given  by 

°s  “ 1 ” cl®2  " C2S  • (13) 


with 


and 


c1  - 0.001982 


c2  * 0.00036 

Equations  (1)  - (4)  are  solved  at  sigma  levels  defined  by 
setting  s ■ 1,  2,  3,..*  20  in  (13)  while  the  eddy  diffu- 
sion coefficient  is  calculated  at  levels  defined  by 
s * 3/2,  5/2,  7/2,..,  39/2.  This  choice  of  parameters  places 
the  first  grid  point  at  approximately  20  m above  the  sea 
surface  and  the  last  grid  point  at  about  tropopause  level. 

The  numerical  scheme  is  forward  in  time  with  the 
eddy  flux  terms  in  (1)  - (4)  treated  implicitly.  The  implicit 
treatment  allows  a much  longer  time  step  (typically  5-10  min) 
than  would  otherwise  be  possible. 
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Section  3 

DESCRIPTION  OF  THE  EXPERIMENTS 

Values  for  the  parameters  listed  in  Table  1 (which 
are  considered  to  be  averaged  horizontally  across  the  model 
domain  and  in  time  over  a period  of  about  one  month)  are 
imposed,  the  sea  surface  temperature  is  specified,  and 
(1)  - (4)  are  integrated  to  a steady  state  from  arbitrary 
initial  conditions.  The  resulting  solutions  represent  the 
mean  state  of  the  atmosphere  averaged  in  time,  and  space 
as  above.  Then,  with  the  same  values  for  the  imposed  para- 
maters  in  Table  1.  the  model  is  integrated  to  other  steady 
states  corresponding  to  various  sea,  surface  temperature 
perturbations,  with  the  resulting  differences  in  the 
thermodynamic  characteristics  of  the  atmosphere  being  the 
quantities  of  interest. 

Three  experiments  are  performed  with  the  sea  surface 
temperature  specified  as  below: 


Experiment 

Sea  Surface  Temperature 

1 

302.6 

2 

301 .5 

3 

301.0 

These  temperatures  are  in  the  range  characteristic  of  the 
summertime  tropical  ocean. 

Note  that  these  experiments  are  of  a very  preliminary 
nature  and  serve  primarily  to  illustrate  the  basic  capabili- 
ties of  this  very  simple  model.  In  future  studies,  specific 
seasons  and  geographical  locations  will  be  investigated.  In 
these  cases,  the  parameters  in  Table  1 will  have  to  be  chosen 
very  carefully  from  climatological  data  rather  than  specified 
in  a simple  fashion  as  is  done  here.  In  addition,  the  present 
model  lacks  a parameterization  for  clouds  and  rainfall.  These 
effects  are  obviously  fundamental  in  modeling  the  atmosphere 
in  certain  regions  and  during  certain  periods  and  must  there- 
fore be  Incorporated  into  later  versions  of  the  model. 
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Table  1 

Imposed  Parameters  For  All  Three  Experiments 


Parameter 

E-W  component  of  geostrophic  wind 

Symbol 

Ug 

Value 

10  m s"1 

N-S  component  of  geostrophic  wind 

Vg 

i 

0 

specific  humidity  at  distance  d 
upstream 

qus 

see  Fig.  1 

1 

length  scale  for  lateral  moisture 
advection 

d 

300  km 

surface  pressure 

P* 

1000  mb 

vertical  motion  in  o coordinate 
system 

e 

o 

see  Fig.  1 

specific  humidity  at  model  top 

qt 

0.09  g kg  _1 

temperature  at  model  top 

Tt 

220°K 

E-W  component  of  wind  at  model  top 

ut 

10  m s"1 

N-S  component  of  wind  at  model  top 

vt 

0 

latitude 

♦ 

25° 

declination  angle  of  sun 

6 

20° 

Qus  (9  kg”1)  0 p*  (mb  hr"1) 

Figure  1.  Prescribed  profiles  of  the  upstream  humidity  q and  the 
vertical  motion  0 p*  in  the  lowest  1500  m,  forus 
•11  three  experiments. 
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section 4 
RESULTS 


Figure  2 shows  vertical  profilesof  potential  tem- 
perature and  specific  humidity  in  the  lowest  1500  m of  the 
atmosphere  for  the  steady  state  attained  in  experiment  1. 

Note  the  presence  of  (a)  the  surface  layer  of  thickness 
20  m in  which  potential  temperature  and  specific  humidity 
decrease  rapidly  with  height  and  (b)  the  Ekman  layer  of 
thickness  400  m in  which  potential  temperature  and  specific 
humidity  vary  almost  linearly  and  show  only  a slight  decrease 
with  height.  The  Ekman  layer  is  topped  by  a stable  layer 
in  which  the  specific  humidity  decreases,  first  rapidly 
and  then  gradually,  with  height.  These  results  represent 
a reasonable  thermodynamic  state  for  the  cloud-free  marine 
atmosphere. 

Figure  3 shows  vertical  profiles  of  u and  v for  the 
steady  state  of  Experiment  1.  Note  the  large  value  of 
wind  shear  in  the  surface  layer  and  the  clockwise  rotation 
of  the  wind  vector  with  height  in  the  Ekman  layer,  both  of 
which  are  expected  results. 

The  vertical  profiles  of  u and  v for  Experiments  2 
and  3 are  essentially  the  same  as  in  Experiment  1.  However, 
the  lowest  500  m in  Experiment  2 is,  on  the  average,  0.4°C 
cooler  and  0.3  g kg*1  dryer  than  Experiment  1,  while  in 
Experiment  3 this  layer  is  0.8°C  cooler  and  0.6  g kg*1  dryer 
than  Experiment  1.  Hence,  as  expected,  lower  sea  surface 
temperatures  produce  a cooler  and  dryer  marine  boundary  layer. 
Furthermore,  for  the  range  of  temperature  perturbations 
considered,  the  decrease  in  the  average  temperature  and 
moisture  in  the  lower  atmosphere  is  a linear  function  of. the 
sea  surface  temperature. 
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Table  2 gives  values  for  the  surface  fluxes  of  latent 
heat,  sensible  heat,  long-wave  radiation,  and  short-wave 
radiation  for  all  three  experiments.  The  total  surface  heat 

balance,  E , for  Experiment  1 indicates  a net  heating  rate 

0 -2  -1 
for  the  upper  ocean  of  102  cal  cm  day  , The  net  heating 

-2  —1 

rates  for  Experiments  2 and  3 are  118  and  133  cal  cm  day 
respectively.  Hence,  the  model  predicts  a linear  variation 
of  Eq  with  the  sea  surface  temperature  for  perturbations  of 
1°C  or  less.  Furthermore,  we  find  dE  /dT  = 31  cal  cm-^ 
day-  °C-  which  is  about  half  the  value  estimated  by  Bathen 
et  al.  (1976)  for  a proposed  OTPP  site  at  Keahole  Point, 
Hawaii,  and  about  one-third  the  value  estimated  by  Piacsek 
et  al.  (1976)  for  a proposed  site  off  the  southeast  coast 
of  Puerto  Rico. 

Both  Piacsek  et . al.  (1976)  and  Bathen  et  al.  (1976) 
used  bulk  aerodynamic  formulations  and  local  meteorological 
data  to  estimate  the  surface  fluxes,  and  were  forced  to  make 
certain  arbitrary  assumptions  about  the  atmosphere's  adjust- 
ment to  the  sea  surface  temperature  perturbations.  By 
modeling  the  atmosphere  and  explicitly  predicting  its 
adjustment  to  the  sea  surface  temperature  perturbation,  we 
feel  a more  reliable  determination  of  dEQ/dTo  can  be  obtained. 


Table  2 

Variation  of  the  Surface  Heat 
Flux  with  the  Sea  Surface  Temperature 


r 

Experiment 

1 

2 

3 

Sea  surface  temperature 

302. 0°K 

301. 5°K 

301. 0°K 

Net  solar  radiation 

920 

920 

920 

Net  infrared  radiation 

-449 

-446 

-444 

r 

Latent  heat  flux 

-318 

-307 

-296 

Sensible  heat  flux 

-51 

- 49 

-47 

Net  surface  heat  flux 

102 

118 

133 

c 

(Units  for  fluxes  are  cal 

-2  . -1 
cm  day 

Positive 

flux 

indicates  flow  of  heat  from  air  to  water.  ) 
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Section  5 


SUMMARY  AND  CONCLUSIONS 

A one-dimensional  numerical  model  has  been  con- 
structed and  used  to  investigate  the  effect  of  sea  surface 
. temperature  perturbations  generated  by  Ocean  Thermal  Power 
Plants  on  the  thermodynamic  structure  of  the  cloud-free  marine 
atmosphere.  Conservation  equations  for  momentum,  moisture, 
and  heat  are  solved  on  a vertically  stretched  grid  extending 
from  the  sea  surface  to  tropopause  level.  Long  and  short 
wave  radiation,  horizontal  and  vertical  advection  of  heat 
and  moisture  by  the  large  scale  flow  and  vertical  eddy  fluxes 
of  momentum,  moisture,  and  heat  are  included.  The  vertical 
eddy  fluxes  are  determined  from  Monin-Obukhov  similarity 
theory  for  the  surface  layer,  while  in  the  remainder  of  the 
planetary  boundary  layer,  a semi-empirical  mixing  length 
formulation  is  used. 

Appropriately  averaged  values  of  the  geostrophic 
wind  and  the  large  scale  vertical  motion  field  are  imposed 
along  with  boundary  conditions  on  the  prognostic  variables, 
and  the  model  is  integrated  to  a steady  state.  Then,  for 
small  changes  in  the  sea  surface  temperature  but  with  all 
other  imposed  parameters  remaining  the  same,  other  steady 
states  are  determined  with  the  resulting  differences  in  the 
thermodynamic  characteristics  of  the  atmosphere  being 
the  quantities  of  interest. 

It  is  found  that  for  sea  surface  temperature 
perturbations  of  1°C  or  less,  the  average  temperature  and 
moisture  in  the  marine  boundary  layer  vary  almost  linearly 
with  the  sea  surface  temperature  with  slopes  of  0.8  °C/°C 
and  0.6  g kg_1/°C  respectively.  The  total  heat  flux  at  • 
the  sea  surface  is  also  found  to  vary  linearly  with  the 
sea  surface  temperature  in  the  range  with  a slope  of  31  cal 
cm”*5  day”1/°C. 
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Although  these  preliminary  results  indicate  that 
OTPP  operation  will  induce  only  a very  minor  change  in  the 
thermodynamic  structure  of  the  marine  atmosphere,  it  must 
be  reemphasized  that  the  effect  of  cumulus  clouds  is  ignored 
in  the  model.  Hence,  the  predicted  changes  in  the  atmos- 
- phere's  state  for  various  sea  surface  temperature  perturba- 
tions apply  only  to  cloud-free  conditions  unless  the 
assumption  is  made  that  the  characteristics  of  the  overlying 
cumulus  cloud  ensemble  are  unaffected  by  the  OTPP-generated 
sea  surface  temperature  anomaly.  However,  this  is  probably 
not  a good  assumption.  For  example,  the  observations  re- 
ported by  Malkus  (1957)  suggest  that  trade  wind  cumuli  respond 
to  sea  surface  temperature  anomalies  as  small  as  0.1  - 0.3°C. 
Therefore  changes  induced  in  the  thermodynamic  structure  of 
the  sub-cloud  layer  by  OTPP  operation  are  likely  to  affect 
the  cumulus  clouds  which,  in  turn,  affect  the  thermodynamics 
of  the  sub-cloud  layer  themselves  (Ogura  and  Cho , 1973; 

Esbensen , 1975).  Thus,  feedback  effects  are  possible  that 
might  amplify  the  surface  flux  perturbations  resulting  from 
OTPP  operation. 

Furthermore,  because  of  the  tendency  for  equivalent 
potential  temperature  to  be  conserved  in  moist  convection 
and  the  exponential  dependence  of  saturation  specific  humidity 
on  temperature,  small  changes  in  the  temperature  at  cloud 
base  yield  much  larger  changes  in  the  temperature  at  cloud 
top.  Therefore,  to  achieve  a better  prediction  of  the  atmos- 
phere's adjustment  to  OTPP-generated  sea  surface  temperature 
anomalies,  the  model  should  be  extended  to  include  a 
parameterization  of  the  thermodynamic  effect  of  cumulus  clouds. 
The  cumulus  parameterization  scheme  of  Arakawa  and  Schubert 
(1974)  provides  the  only  closed  theory  for  the  mutual  inter- 
action of  a cumulus  cloud  ensemble  with  the  large-scale 
environment  and  is  now  being  incorporated  into  the  model. 
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Appendix 
List  of  Symbols 


B 

c_ 


C1  ’ c2 


T*  F* 
*L,fS 


albedo  of  sea  surface 

Stef an-Boltzmann  constant 

specific  heat  at  constant  pressure  for  air 

stretched-grid  parameters 

horizontal  length  scale  for  calculation  of 
lateral  moisture  advection 

net  heat  flux  at  sea  surface  (positive  downward) 

net  downward  flux  of  long-wave  radiation, 
short-wave  radiation 


f Coriolis  parameter 

g acceleration  of  gravity 

H ,H  vertical  eddy  flux  of  sensible  heat  in  surface 

layer,  Ekman  layer 

K vertical  diffusion  coefficient  for  momentum, 

moisture,  and  heat 

£ mixing  length 

L Monin-Obukhov  length 

N Brunt-Vaisala  frequency 

P»P*.P  , Prp  pressure,  surface,  pressure,  reference  pressure, 

r pressure  at  model  top 

Q , Q vertical  eddy  flux  of  moisture  in  surface  layer, 

Ekman  layer 

q,q  specific  humidity,  specific  humidity  at  sea 

° surface 

q specific  humidity  at  distance  d upstream  from 

us  the  center  of  the  model  domain 


A-l 


c 


Tv,Tr 


u,u„,u„ 
g go 


v’  Vg,Vgo 


Vwt 


gas  constant  for  dry  air 
solar  constant 

temperature,  sea  surface  temperature 
virtual  temperature,  reference  temperature 
wind  speed,  wind  speed  at  sea  surface 

geostrophic  wind  at  top  of  Ekman  layer,  (u  2 + 

x ' 

surface  friction  velocity  (t  /p  )“ 

0 ft 

E-W  component  of  wind,  E-W  component  of  geo- 
strophic wind,  E-W  component  of  geostrophic 
wind  at  sea  surface 

N-S  component  of  wind,  N-S  component  of 
geostrophic  wind,  N-S  component  of  geostrophic 
wind  at  sea  surface 

precipitable  water  vapor 

precipitable  water  vapor  at  sea  surface, 
model  top 

height  above  sea  surface 

height  at  which  Richardson  number  equals  1 
surface  roughness  length 


pa-  pw 
Vvw 


height  of  first  grid  point  above  sea  surface 
von  Harman’ s constant 
density  of  air 

reference  density  for  air,  water 
molecular  viscosity  of  air,  water 
latitude 

non-dimensional  gradients  of  wind  speed, 
temperature,  and  humidity  in  the  surface  layer 

wind  stress  at  sea  surface 


Q* 


1 


VTy 


a 


© 


e 


C! 


vertical  flux  of  E-W  momentum,  N-S  momentum 
(i.e.,  components  of  turbulent  shear  stress) 

vertical  coordinate,  p/p* 

vertical  motion  in  o coordinate  system,  d /dt 

-R/c_ 

potential  temperature,  To  * 

infrared  emissivity  of  sea  surface 

derivative  of  Infrared  emissivity  with  respect 
to  precipitable  water  vapor 

solar  zenith  angle 


A-3 
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